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ABSTRACT 

Context. Blue supergiant stars are known to display photometric and spectroscopic variability that is suggested to be linked to stellar 
pulsations. Pulsational activity in massive stars strongly depends on the star’s evolutionary stage and is assumed to be connected with 
mass-loss episodes, the appearance of macroturbulent line broadening, and the formation of clumps in the wind. 

Aims. To investigate a possible interplay between pulsations and mass-loss, we carried out an observational campaign of the supergiant 
55 Cyg over a period of five years to search for photospheric activity and cyclic mass-loss variability in the stellar wind. 

Methods. We modeled the H, He i. Si ii, and Si iii lines using the nonlocal thermal equilibrium atmosphere code FASTWIND and de¬ 
rived the photospheric and wind parameters. In addition, we searched for variability in the intensity and radial velocity of photospheric 
lines and performed a moment analysis of the line profiles to derive frequencies and amplitudes of the variations. 

Results. The Ha line varies with time in both intensity and shape, displaying various types of profiles: P Cygni, pure emission, almost 
complete absence, and double or multiple peaked. The star undergoes episodes of variable mass-loss rates that change by a factor of 
1.7-2 on different timescales. We also observe changes in the ionization rate of Sin and determine a multiperiodic oscillation in the 
He I absorption lines, with periods ranging from a few hours to 22.5 days. 

Conclusions. We interpret the photospheric line variations in terms of oscillations in p-, g-, and strange modes. We suggest that these 
pulsations can lead to phases of enhanced mass loss. Furthermore, they can mislead the determination of the stellar rotation. We 
classify the star as a post-red supergiant, belonging to the group of a Cyg variables. 

Key words. Stars: early-type - Stars: supergiants - Stars: winds, outflows - Stars: mass-loss - Stars: activity - Stars: individual: 
55 Cyg 


1. Introduction 

Blue supergiants (BSGs) are evolved, luminous objects that 
often display strong photometric and spectroscopic variability 
(e.g., Rosendhal 1973; Kaufer et al. 1997; Waelkens et al. 1998; 
Mathias et al. 2001; Kaufer et al. 2006; Lefever et al. 2007; Clark 
et al. 2010) Long-term space-based photometry has linked this 
variability to stellar pulsations, establishing a new instability do¬ 
main in the Hertzsprung-Russell (HR) diagram region populated 
by BSGs (Saio et al. 2006) The presence of pulsations in BSGs 
requires that excited modes are reflected either at an intermediate 

* Based on observations taken with the Perek 2m telescope at Ondfe¬ 
jov Observatory, Czech Republic, and the Poznan Spectroscopic Tele¬ 
scope 2 at the Winer Observatory in Arizona, USA 


convection zone connected to the hydrogen-burning shell (Saio 
et al. 2006; Godart et al. 2009), or at the top of the chemical 
composition gradient region surrounding the radiative He core 
(Daszyhska-Daszkiewicz et al. 2013) 

While BSGs in the vicinity of the main sequence are found 
to pulsate in both pressure (p-modes with periods of hours to 
about one day) and gravity modes (g-modes with periods of 2- 
10 days), objects slightly more evolved to the red are predicted to 
pulsate in pure g-modes (Saio et al. 2006) However, Kraus et al. 
(2012) recently discovered a periodic variability of 1.6 hours in 
the late-type BSG star HD 202 850. Existing pulsation models 
do not predict such short periods; this emphasizes the current 
deficiency in our knowledge of the pulsation activity in BSGs. 
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As massive stars can cross the BSG domain more than once, 
the pulsational activity of these stars can drastically change be¬ 
tween their red- and blueward evolution (Saio et al. 2013) When 
compared with their less-evolved counterparts, BSGs on a blue 
loop, or blueward evolution, tend to undergo significantly more 
pulsations, even including radial strange-mode pulsations (with 
periods of 10-100 or more days), as observed in a Cygni vari¬ 
ables. 

Additionally, the photospheric line profiles of BSGs typ¬ 
ically contain some significant extra broadening, so-called 
macroturbulence (Abt 1958; Conti & Ebbets 1977), which is 
often on the same order as rotational broadening (Ryans et al. 
2002; Simon-Dfaz & Herrero 2007; Markova & Puls 2008) Re¬ 
cent evidence links macroturbulence and line profile variability 
to stellar pulsational activity (Aerts et al. 2009; Simon-Diaz et al. 
2010a,b) providing additional proof for the presence of pulsa¬ 
tions in BSGs. 

BSGs lose mass via their line-driven winds, but the rates at 
which the material is lost are still highly debated. Diverse spec¬ 
tral ranges trace distinct wind regions, and comparison of mass- 
loss rates obtained from various methods typically results in val¬ 
ues that can differ by a factor 2-10. In addition, the observed 
complex wind structures seen in the Ha variability of BSGs (e.g., 
Kaufer et al. 1996; Markova et al. 2008) for instance, imply that 
these objects cannot have smooth spherically symmetric winds. 
Spectropolarimetric observations revealed no evidence of mag¬ 
netic fields (Shultz et al. 2014) Instead, the winds of BSGs are 
found to be clumped (Prinja & Massa 2010) Micro- and macro- 
clumping severely influence both mass-loss determination and 
line profile shapes of not only resonance lines, but also Ha (e.g., 
Sundqvist et al. 2010, 2011, 2014; Surlan et al. 2012, 2013) 
The onset of clumping in the winds of massive stars could be 
due to line-driven instability (see, e.g., Feldmeier 1998) Pertur¬ 
bations leading to this instability might be initiated by pulsa¬ 
tional activity, at least in the vicinity of the stellar photosphere. 
Moreover, strange mode pulsations have been suggested to cause 
time-variable mass loss in luminous, evolved massive stars (see 
Glatzel et al. 1999; Aerts et al. 2009; Puls et al. 2011) The first 
observational evidence of these modes was found in the massive 
BSG HD 50064 (Aerts et al. 2010b) 

The ability of BSGs to maintain stable pulsations opens a 
completely new perspective in studying these stars. Based on 
established methods from asteroseismology, it will be possible 
to investigate not only the stellar atmospheres from which the 
mass-loss and possible onset of wind clumping are initiated, but 
also the deep interior of the stars, revealing important physi¬ 
cal properties, such as internal structure, rotation, and mixing. 
Knowledge of these parameters is of vital importance for our 
understanding of the post-main sequence evolution of massive 
stars, as these are key input parameters to modern stellar evo¬ 
lution calculations. This emphasizes the need for further inves¬ 
tigations, both theoretically and observationally, to improve our 
understanding of pulsational activity in BSGs in different evolu¬ 
tionary stages. 

To investigate, in particular, the possible interplay between 
pulsations and mass-loss, we initiated an observational cam¬ 
paign to search for pulsational activity and cyclic mass-loss vari¬ 
ability in stellar winds for a sample of bright BSGs. Here we 
report on our results for the object 55 Cyg. 

2. The star 

55 Cyg (HD 198478, MWC353, HR 7977, HIP 102 724, a = 
20''48”56?291 and 6 = -h 46° 06' 50788; 2000) is a bright star 


located in the Cyg OB7 association (Humphreys 1978) that was 
classified as B3Ia by Morgan & Roman (1950) Its hipparcos 
parallax of 1.40 + 0.17 mas (van Leeuwen 2007) places the star 
at a distance of 714 ± 86 pc. The interstellar visual extinction to¬ 
ward 55 Cyg was found to be 1.62 mag (Barlow & Cohen 1977; 
Humphreys 1978) 

Studies of its photospheric lines revealed two important 
facts: different elements typically display different radial veloc¬ 
ities (e.g., Hutchings 1970) while the same lines show strong 
variations with time. For instance, Underhill (1960) described 
radial velocity variations of 25 kms with no evidence that the 
changes found were periodic. Eater, Granes (1975) collected 34 
spectra, distributed over 15 consecutive nights, and found that 
the radial velocity curves oscillate with a period of 4-5 days, in 
addition to photospheric motions on timescales about three times 
longer. 

Photometric observations are found to display microvari¬ 
ability, but clear periodicities are difficult to establish. While 
Rufener & Bartholdi (1982) reported an 18-day variability in 
their V-band magnitudes, Koen & Eyer (2002) found a period¬ 
icity of 4.88 days using hipparcos data. On the other hand, no 
indication for clear periodicities were seen by Percy & Welch 
(1983) van Genderen (1989) and Lefevre et al. (2009) 

The Ha line in 55 Cyg is also highly variable. While Under¬ 
hill (1960) and Rzaev (2012) remarked that in their spectra the 
Ha line maintained a P Cygni profile, but with highly variable 
strength, Ebbets (1982) noted that during one of his observa¬ 
tions, Ha showed additional faint emission wings, which were 
not visible anymore in the spectrum he took two months later. 
Furthermore, Maharramov (2013) discovered that Ha can even 
completely disappear from the spectrum. 

The stellar rotation velocity, projected on the line of sight 
(v sin i), and the terminal wind velocity (Uoo) are also highly vari¬ 
able or uncertain. For both, broad ranges are found in the litera¬ 
ture. Values for v sin i range from 0 km s ' (Slettebak & Howard 
1955) to a maximum of 61kms ' (Howarth et al. 1997) The 
latter was obtained from lines in the lUE spectra. Other de¬ 
terminations based on optical spectral lines and different tech¬ 
niques include 42kms“' (Day & Warner 1975), 49 + Ikrns^* 
(Gies & Lambert 1992), 45 + 20kms“‘ (McErlean et al. 1999) 
35 + 9kms“' (Abt et al. 2002) , and 39kms“' (Markova & Puls 
2008) The terminal wind velocity is much harder to determine. 
Based on different lines in the lUE spectra, Prinja et al. (1990) 
obtained 470kms * from the Civ line, while Prinja & Massa 
(2010) derived a higher value of 560 km s ' from Siiv. Both 
values are lower than predictions from line-driven wind theory, 
however, which delivers a value of ~ 690 km s ' (see, e.g., Kr- 
ticka & Kubat 2001) Conversely, Markova & Puls (2008) noted 
that the observed value of 470kms“^ is too high to give satis¬ 
factory fits to their Ha line, which would agree much better with 
a terminal wind velocity of only ~ 200kms '. These authors 
speculated that 55 Cyg might have a time-variable wind veloc¬ 
ity. 

There is also scatter in the derived stellar parameters. Al¬ 
though the determination of fundamental parameters could be bi¬ 
ased by the technique used, different values of are found us¬ 
ing models built on similar approaches. Teff measurements based 
on the Balmer discontinuity give values of 16 450 K (the color 
index method of Gies & Lambert 1992) and 15 380 + 860 K (the 
BCD spectrophotometric method of Zorec et al. 2009) respec¬ 
tively. Using nonlocal thermal equilibrium (non-LTE) plane- 
parallel hydrostatic stellar atmosphere and the Si ionization bal¬ 
ance, the derived Tetr values are 18 000 K (McErlean et al. 1999) 
18 500K (Monteverde et al. 2000) and 17 000K (Jurkic et al. 
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2011) while non-LTE line- and wind-blanketed model atmo¬ 
spheres provide values of 16 500 K (Crowther et al. 2006) and 
17 500 K (Searle et al. 2008; Markova & Puls 2008) In all cases, 
the surface gravity determinations are around log g = 2.2 + 0.1 
(cf. Jurkic et al. 2011) 

As Ha is a fundamental mass-loss diagnostic (e.g., Puls et al. 
1996) it can be expected that model fits to observations ob¬ 
tained at different periods (with correspondingly varied shape 
and line strength) result in different mass-loss rates. And the 
values clearly vary by a substantial amount. While Searle et al. 
(2008) found M - 5x 10“^ Moyr^*, the value of Crowther et al. 
(2006) is a factor of two lower, even though both works used the 
same terminal wind velocity value of 470km s“^ In addition, 
Markova & Puls (2008) obtained a rather large uncertainty in 
their mass-loss values (1.175-4.07) x 10“^ Moyr^', which is due 
to the high uncertainty in the wind terminal velocity (200 km s ’ 
vs 470km s“’). All of these investigations used smooth (i.e., un¬ 
clumped) winds. 

Concerning the possible evolutionary stage of 55 Cyg, Un¬ 
derhill (1969) mentioned that the observed Hy line profile was 
too shallow, which could indicate that 55 Cyg is H-poor, and 
hence a post-red supergiant. A similar conclusion was drawn by 
Lennon et al. (1993) Gies & Lambert (1992) Crowther et al. 
(2006), and Searle et al. (2008) based on the CNO-processed 
material seen in its optical spectra, in particular, its strong N- 
enrichement. Markova & Puls (2008) additionally concluded 
that the star probably is enriched in He. These authors also found 
a rather low spectroscopic mass of ~ 11 Mq compared to the 
star’s high luminosity, which places it on an evolutionary track 
of an initial mass of 25 Mq. Additional support for its evolved 
nature comes from van Genderen (1989) who listed 55 Cyg as 
an a Cygni variable (i.e., a luminous, photometrically variable 
supergiant of spectral type B and A). If this classification is cor¬ 
rect, then 55 Cyg should display numerous pulsation modes, in¬ 
cluding radial strange modes (Saio et al. 2013) 

3. Observations 

We spectroscopically monitored 55 Cyg between 2009 August 
15 and 2013 October 22. We obtained a total of 344 spectra, dis¬ 
tributed over 64 nights, with the Coude spectrograph attached 
to the Perek 2 m telescope at Ondfejov Observatory (Slechta & 
Skoda 2002) Until the end of May 2013, the observations were 
taken with the 830.77 lines mm“’ grating and a SITe 2030 x 800 
CCD. Beginning in June 2013, we used the newly installed Py- 
LoN 2048 X 512BX CCD. With both detectors, a spectral resolu¬ 
tion of R ^ 13 000 in the Ha region was achieved, and the wave¬ 
length coverages were from 6253 A to 6764 A for the old CCD 
and from 6263 A to 6744 A with the new one. For wavelength 
calibration, a comparison spectrum of a Th-Ar lamp was taken 
immediately after each exposure. The stability of the wavelength 
scale was verified by measuring the wavelength centroids of O i 
sky lines. The velocity scale remains stable within 1 kms ’. In¬ 
dividual exposure times range from 400 s to 2700 s, and, for all 
spectra, the achieved signal-to-noise (S/N) ratio was higher than 
300. During the nights of 2013 July 19-22 and 2013 August 10- 
12, we achieved S/N values of more than 500, allowing for high- 
quality time series, covering several hours of continuous obser¬ 
vation. 

We also obtained a series of echelle spectra using the Poz¬ 
nan Spectroscopic Telescope at the Winer Observatory in Ari¬ 
zona, USA. Observations were carried out from 2013 October 
17 to 2013 November 7, and a total of 41 spectra were acquired. 



Fig. 1: Spectrum of 55 Cyg taken in the Ha region with the Perek 
2m telescope at Ondfejov Observatory. 
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Fig. 2: Night-to-night variation in the Ha line in the observing 
period 2013 July 19 to August 12. This period shows all of the 
typical profile shapes: P Cygni, pure emission, almost complete 
absence, and double and multiple peaked. 

covering the wavelength range 3900 to 8677 A, with a spectral 
resolution of = 40 000. Each spectrum has an exposure time 
of 1800 s and a S/N of ~ 120- 180. 

All data were reduced and the heliocentric velocity corrected 
using standard IRAfQtasks. We also observed a rapidly rotating 
star (HR 7880, Regulus, ^Aql) once per night to perform the 
telluric correction. 

A typical Ondfejov spectrum is shown in Fig.[2 In general, 
the overall appearance of the full optical spectrum is very similar 
to the ones presented by Lennon et al. (1992) and Chentsov & 
Sarkisyan (2007) Interesting is the shallow, but broad emission 
component centered on the Ha line, which is noticeable in each 
of our spectra, while it was reported only once in the literature 
(Ebbets 1982) The extent of the emission wings measured from 
our spectra is +1200kms ’. Such broad emission components 
typically originate from incoherent electron scattering (see, e.g., 
McCarthy et al. 1997; Kudritzki et al. 1999) 

The Ha line displays time variations in both intensity and 
shape, showing various types of profiles: P Cygni, pure emis¬ 
sion, almost complete absence, and double or multiple peaked. 

' IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation. 
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As an example, see Fig.|^ Similar variability of the Ha line has 
been observed in other BSGs that belong to the group of a Cygni 
variables (e.g., Kaufer et al. 1996) 

4. Results 

4.1. Line profile modeling 

To quantify and analyze the spectral line variation of 55 Cyg, we 
simultaneously fit observed line profiles of several representa¬ 
tive transitions of H, He and Si atoms with synthetic line models 
computed with the code Fast Analysis of STellar atmospheres 
with WINDS (FASTWIND). This is a non-LTE, spherically sym¬ 
metric model atmosphere code that enables computating contin¬ 
uum and line radiation fluxes and provides a set of parameters 
that describe the photospheric and wind structure (Santolaya- 
Rey et al. 1997; Puls et al. 2005) The code includes line block¬ 
ing and blanketing effects and consistently calculates the tem¬ 
perature stmcture (see details in Kubat et al. 1999) The wind 
velocity distribution is computed using the so-called yS-law. 

For the fundamental stellar parameters (Tetr, log g) we 
mainly used the Arizona spectra. They cover the Hy and Hd lines 
(known as excellent surface gravity sensors) and the He i A 4471 
triplet line (considered a good gravity indicator, with some sen¬ 
sitivity to Teff as well, Lefever et al. 2007) In addition, we mod¬ 
eled the singlet transitions of He id 6678 A and He id 4713 A, 
and the lines of Siiidd4128,4130A, and Siiiid4552A. These 
Si lines are also excellent temperature indicators (Markova & 
Puls 2008) Unfortunately, the Ondfejov spectra cover a more 
restricted wavelength interval, providing only the He id 6678 A 
line. Although the Ondfejov spectra have a lower resolution, they 
have significantly higher S/N ratios than the Arizona spectra. 

Another input model parameter is the stellar radius. This 
value should be derived in an independent way, such as the de- 
reddened absolute magnitude or the spectral energy distribution. 
Different authors used these methods to derived the radius of 
55 Cyg and obtained values ranging from 38 Rq (Markova & 
Puls 2008) to 83 Rq (Gies & Lambert 1992) Because of the 
large discrepancy in the literature values, we employed other 
methods. Using the BCD parameters, based on a reddening- 
independent method (Zorec et al. 2009) we obtain a bolometric 
magnitude of-8.5+0.3 mag, which leads to log L/Lq = 5.3+0.1 
and R - 61 ± 8Ro- This bolometric magnitude agrees with 
the determination reported by Barlow & Cohen (1977), who 
estimated My = -6.8 mag from the lO/rm excess emission. 
With the bolometric correction of 1.7 (Flower 1996) we derive 
Mboi = -8.6 mag, in agreement with the BCD method. The ra¬ 
dius obtained from these two methods also agrees with the value 
of 61 Rq, derived by Pasinetti Fracassini et al. (2001) based on 
measurements of the stellar angular diameter from infrared pho¬ 
tometry by Blackwell & Shallis (1977) 

The wind parameters were derived by modeling the Ha line. 
To compute the synthetic line profiles for 55 Cyg, we started 
from a set of initial values for both photospheric and wind pa¬ 
rameters, taken from the literature (see Sect. [^. However, these 
values did not deliver good fits, so that we extended the range 
of parameters, looking for the best fit in a “by-eye” procedure. 
For this, we treated all model parameters (Teff, log g, /3, M, Uoc, 
rnaa-o, and Urot) as free values and computed a large grid of syn¬ 
thetic line profiles. The parameters Umicro and R* were initially 
fixed at 10 km s ' and 61 R© , respectively, and then were slightly 
varied to obtain improved fits. We selected the best fit to the ob¬ 
servations and verified the consistency of some of the individual 
parameters, which are not independent variables. For example. 


the stellar radius and the surface gravity should obey the rela¬ 
tion gRl ~ constant. To fit the terminal velocity, we preferen¬ 
tially considered the values that reproduce the blueshifted ab¬ 
sorption component of the Ha P Cygni profile. Otherwise, for 
profiles with pure emission, we derived the terminal velocity 
based mainly on the Ha line width. 

Using the FASTWIND code, we explored the sensitivity of 
the model parameters that characterize both the photosphere and 
the wind. However, we should keep in mind that there are two 
sources of errors. One results from the modeling, in which the 
involved uncertainties result from the statisitcal standard devia¬ 
tion and are obtained using the best-fit model. The second refers 
to the inaccuracy of the input parameters, that is, the radius, re¬ 
sulting from the uncertainties in the stellar distance and hence 
luminosity. But when the distance is fixed, any deviation in the 
parameters larger than the statistical error indicates possible vari¬ 
ability in the physical parameters. 

To determine the statistical uncertainty of the parameters, we 
used the best-fit model and tested different values of each param¬ 
eter, one at a time, keeping all other parameters fixed. We have 
the advantage of the many high-resolution and high S/N ratio 
spectra obtained in Arizona, taken over a sequence of consecu¬ 
tive days. This enabled us to distinguish errors on the mean val¬ 
ues of the studied parameters from real changes of the physical 
conditions. 

From fits to the Sin and Siiii lines we observe uncertain¬ 
ties of 300 K - 500 K in Teff- We also And that the Ha emis¬ 
sion is sensitive to changes in Tetr of ~ 300 K. From the Hy and 
Hd wings, the statistical error in log g is 0.05 dex, but due to 
large uncertainties in other parameters (e.g., the radius), the er¬ 
ror could be slightly larger, but still below 0.1 dex. This latter is 
what we use in the tables as upper limit. The uncertainties in the 
rotational, micro- and macroturbulent velocities depend on the 
modeled line transition. Error bars in (imicro are 2 km s ' for the 
He and Si lines, and 5 km s ' for the H lines, except for Ha, for 
which the error is about 10km s '. We And errors of 5 km s ' 
in iiiot and Umacro for Hy, Hd, He, and Si lines, and of 5 km s ' 
and 10kms ' for Ha. Eor H/l, our error estimates for Uiot and 
iimacro are ~ 20kms“'. To derive the wind parameters, we con¬ 
sidered the shape and intensity of the absorption and emission 
components of the Ha line. We And noticeable line variations 
with changes in M of about 0.1 x lO^’Moyr^' and 10% in Voo 
(although the error in Voo can reach up to 30% if Ha displays 
a pure absorption line). These adopted uncertainties are over¬ 
estimations because they are larger than or equal to the stan¬ 
dard statistical deviation computed with the model parameters of 
the Arizona data, which correspond to a term of 21 days show¬ 
ing small or moderate spectral variations. Unfortunately, even if 
we can check the sensitivity of M by varying the parameters of 
the model, we cannot provide an accurate value of this quantity 
because this depends on the uncertainty of the stellar distance. 
To evaluate the error in relation to the accurate value, we per¬ 
formed an error propagation using the optical depth parameter 
Q - Mf(ucaR*)^'^, which results from the model calculations. 
To compute Q, we selected the observations of Arizona (listed 
in Table[2 log 2 = -12.975 + 0.066) because they were taken 
during many consecutive days and the line proflles did not show 
huge variations. Assuming the scatter in the bolometric magni¬ 
tude AMboi - +0.3, that is, AlogL/Lo ^ 0.12, and the error in 
Teff ^ 500 K, we obtain AlogR»/Ro ^ 0.04. With our uncer¬ 
tainty of 10% in Uoa and the statistical error A log Q ^ +0.07, the 
uncertainty in AM/M ^ 0.3. 

The best-flt values (and corresponding statistical deviations) 
for the photospheric and wind parameters derived for each ob- 
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Fig. 5; Same as Fig.|^ but for data taken in 2011. 
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mined from the Arizona spectra. Both results predict marginally 
higher log g values than those previously reported in the litera¬ 
ture (log g — 2.2 + 0.1), which are typically based on only one 
or two measurements, but the agreement within the error bars is 
fairly good. 

From the 15 Arizona spectra we obtain an average radius of 
57 Rq with a statistical dispersion of +1 This radius is close 
to the initial input value of 61 We noticed that changes in A/?, 
by more than 1 Rq result in appreciable changes in the computed 
prohles. Moreover, several of the 43 Ha profiles in the Ondfejov 
spectra (spread over four years) could not be modeled with a 
fixed radius of 57 + 1 TJq. Instead, best-fit models were obtained 
for values ranging from 52 R^ to 65 Rq. This suggests that the 
star either has a variable wind opacity or undergoes real radial 
changes, for example, via radial pulsations. 

Using the mean value for R* and the surface gravity cor¬ 
rected for rotation, log gftrue = 2.44 dex, we determine a stellar 
mass of 34 ± 4Mo. The mean stellar luminosity is found to be 
logL/Lo = (5.57 + 0.03), which agrees with values derived from 
observations. 


servation are given in Tables [T] (Arizona data) and (Ondfejov 
data). We adopt solar abundances for all elements, but discuss 
the validity, in particular with respect to the He abundance, in 
Sect. l4T4l 

Figure shows a sequence of model fits to simultaneous ob¬ 
servations of Ha, H/3, Hy, Hd, He 1 4 6678 A, He i T 4471 A, and 
He 1 4 4713 A line spectra taken in Arizona during 15 succes¬ 
sive nights in 2013 (see also Table [Tjl . The fits to the lines of 
Si 11444128,4130 A, and Si iii44552 A are shown in Fig.|^ We 
obtained very good fits for the whole sample of lines, with the 
exception of the blue wing of He 1 4 4471 A, where we observe a 
tiny discrepancy between the model and the spectrum. This dis¬ 
crepancy corresponds to the presence of a synthetic weak Hei 
forbidden component at 44470 A, which is not observed in any 
of our spectra. Ha shows variations in strength in both the ab¬ 
sorption and emission components of the P Cygni profile. In¬ 
terestingly, the strength and width of the Si ii lines shows night- 
to-night variations, while the Si ni lines remain almost constant. 
The obtained log (EW(Siii4 4128)/EW (Si iii4 4553)) ranges 
between -0.3 and -0.5. This variation can be attributed to a small 
change in the ionization rate of Si ii. 

Fits to the time-scattered observations taken in Ondfejov are 
displayed in Figs. m As these spectra contain only two lines 
that can be modeled, that is, Ha and He 1 4 6678 A, the scatter 
in the obtained values is slightly larger. In general, the models 
match the Ha line profile very well, with the exception of those 
observations displaying complex features in absorption and/or 
double or multiple peaked emissions. We note an interesting 
event of enhanced mass-loss, observed on 2011 September 25, 
when the Ha line shows an unusual P-Cygni profile with a deep 
absorption component and a small emission feature (see Fig.|^. 

4.1.1. Stellar parameters 

On average, we represent the photosphere of the star with Teff = 
18 800 K, with variations ranging from 18 570 K to 19 100 K (see 
Tables and 1^. We find a mean value for log g of 2.43 with a 
scatter of +0.14 dex, based on the fits of the Hy line profile. The 
same trend is obtained from the Ondfejov spectra based on the 
analysis of the Ha line alone. Although these latter log g values 
can be less accurate, they range from 2.3 to 2.5 dex, but their 
average is the same (2.43 + 0.1 dex) as the mean value deter¬ 


4.1.2. Wind variability 

To gain insight into the variability of the stellar wind of 55 Cyg, 
we modeled the Ha line profile over the observing period of 
about five years. 

We fit the observations assuming a smooth (unclumped) 
wind structure represented by a /^-velocity law, with [i - 2.0, 
except in rare cases, where the value of fi was slightly adjusted 
(1.6 and 2.2). On average, the derived M from the Ondfejov data 
is ~ 2.5 X lO'^Moyr^*, with fluctuations of a factor of three 
between the minimum value of 1.5 x 10“’ Mq yr“’ and the max¬ 
imum value of 4.6 x 10“’ Mq yr“'. Our results agree on average 
with the homogeneous wind parameters reported by Crowther 
et al. (2006): P - 2 and M - 2.3 x 10“’Moyr“*, and our val¬ 
ues for the mass-loss fall between the extreme ones obtained 
for unclumped winds in the literature (see Sect.|^. In addition, 
the variation in the terminal velocity (ranging typically from 
230kms“^ to 350kms“') presents a positive correlation with 
the mass-loss rate. There are two exceptions with substantially 
higher values of 600kms“' and 700kms“\ corresponding to 
the two occasions of highest mass loss of 4.6 x 10“’ Mq yr“' and 
4.0 X 10“’Moyr“*. The corresponding observations were taken 
on 2011 September 25 and 2012 February 26, when a mass erup¬ 
tive event probably occurred. The terminal velocities found are 
higher than the 200kms“' reported by Markova & Puls (2008) 
and lower than the 470 km s“* adopted by Crowther et al. (2006) 

The average mass-loss rate obtained (Ondfejov observations) 
is consistent with the average M value estimated from the high- 
resolution spectra taken in Arizona during October and Novem¬ 
ber, 2013. Furthermore, the Arizona spectra show changes in the 
intensity of the Ha line (as shown in Fig.[^ and variations in the 
mass-loss rate by a factor of 1.9. 

4.1.3. Turbulent and projected rotation velocities 

In addition to the stellar and wind parameters discussed so far, 
three more parameters were used to match the observed line pro¬ 
hles. These are microturbulence, macroturbulence, and projected 
rotation velocity. Microturbulent velocities refer to atmospheric 
motions on scales below the photon mean free path, and there is 
evidence that these motions could result from subsurface convec¬ 
tion caused by the iron opacity peak (Cantiello et al. 2009) The 
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Fig. 3: Sequence of fits to daily line variations of H and He lines in 2013. Fits to the different line transitions are achieved using 

(for details see Sect.|4. 1.3||. 
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Fig. 3: Continued. 
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6550 6570 6550 6570 6550 6570 6550 6570 6550 6570 

Wavelength (A) Wavelength (A) Wavelength (A) Wavelength (A) Wavelengh (A) 


Fig. 4: Ha line profiles taken in 2009 and 2010 at Ondfejov Observatory together with their corresponding model fits (for details 
see Table|^. Lines with highly complex profiles or roughly zero EW have not been fitted. 


values found for BSGs range between 5 km s ' and 15 km s ' 
with a tendency to decrease with decreasing temperature (e.g., 
Markova & Puls 2008) Fits to the helium absorption lines give 
best results for 10 + 2kms ', except for He id 6678, which re¬ 
quires a higher value (17 + 2 km s“'). Similar values are observed 
in Sin and Siiii lines, with 12 + 2kms“' and 7 + 2kms“', re¬ 
spectively. The line profiles of Ha also require a higher micro- 
turbulent velocity to achieve satisfactory fits. Here, a value of 
~ 50 + 5kms"' is needed, which can be explained by the fact 
that Ha is formed in the lower wind where supersonic turbulent 
velocities can occur. Problems in fitting Ha line profiles are not 
new, and, for many BSGs, values as high as 50 km s ' (or tur¬ 
bulent velocities increasing linearly from the photospheric value 
to about 50 km s ') were reported (e.g., Crowther et al. 2006; 
Searle et al. 2008) A possible interpretation of the greater mi- 
croturbulent velocity value for He i A 6678 might be that the line¬ 
forming region extends into the base of the wind where it is af¬ 


fected by the wind turbulence. We return to this point in Sect. 5.2 


After including the microturbulence, the synthetic line pro¬ 
files must be broadened by the combined effect of stellar rotation 
and macroturbulence. Considering that v sin i, as a global param¬ 
eter, should be the same for all lines, the variability in the line 
profiles must be assigned to macroturbulent broadening. For H 
and He I lines, the best fit was achieved with v sin / = 55 + 5 


km s“' (see also Sect. 5.2 1 , but theoretical and observed Si ii and 
Sim lines match better if iisin/ is taken around 40 + 5 kms^' 
or 45 + 5 km s '. We also note that the macroturbulence can 
be quite different in various lines. For instance, from modeling 
the Arizona data, we found that fmacro is relatively small in the 
lines of He IT 4471 (6-7 km s-'), He 1 4 4713 (11-12 km s-'), and 
H5 (14-15 kms“'), and approximately constant throughout the 
observing period. In Hy, the changes are larger (14-30 km s '). 
The situation changes for and He id 6678, for which much 
higher and strongly variable values are found (25-90 km s ' and 
20-50 km s ', respectively). Those obtained for Ha vary in the 
range 20-60 km s’ (see Table [^. The Sin and Sim lines also 
show a large dispersion in Umacro, they range from 10 kms^' to 
30 kms^', and 30 kms“ 
tanties of + 5 kms^’. 


to 50 km s ’, respectively, with uncer- 


The data from Ondfejov provide macroturbulent velocities 
only for Ha and He i A 6678 A. Over the entire observing period, 
we obtain values in the range 30-60 km s ’ for He id6678 A, 
while for Ha, occasionally values of more than 100 km s ' are 
found (see Table |^. In most spectra, we find that for both lines 
the contribution of macroturbulence to the line broadening is 
comparable to or higher than that from u sin i. Such high values 
are commonly found in BSGs (e.g., Ryans et al. 2002; Simon- 
Dfaz & Herrero 2007; Markova & Puls 2008; Markova et al. 
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Fig. 6: Same as Fig.|^ but for data taken in 2012 (top) and 2013 (bottom). 


2014; Simon-Dfaz & Herrero 2014) and the values of v sin i and 
Umacro found for 55 Cyg agree well with those of BSGs of compa¬ 
rable mass and effective temperature (see Markova et al. 2014) 


4.1.4. Atmospheric chemical composition 

The Ha line is known to be sensitive to variations in the He con¬ 
tent, in the sense that He-rich wind models produce Ha EWs 
that are a factor of 2-3 lower than those with solar He abundance 
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Table 1: Stellar and wind parameters obtained from line profile fits of observations taken in Arizona in 2013. Listed are the date of 
observation, the Julian date (JD), the derived effective temperature (Teff) and log g, the {3 exponent of the velocity law, the mass-loss 
rate (M), the terminal velocity {Voo), the macroturbulent velocity in Ha (i'macro(Ha)), the stellar radius (J?*), and the luminosity (L). 
Mass-loss rates are for smooth wind models. 


Date 

dd/mm/yr 

JD 

Teff 

[K] 

log g 
[cgs] 

P 

M 

[10-®Moyr-'] 

t^oo 

[kms-'l 

f^macro(Ho^) 

[kms-1] 

R* 

[Re] 

log L/Lo 

17/10/2013 

2456582.72 

19 000+400 

2.50+0.1 

2.0 

0.295+0.01 

310+30 

50+5 

58+1 

5.60+0.03 

18/10/2013 

2456583.70 

19 000+400 

2.50+0.1 

2.0 

0.295+0.01 

310+30 

20+5 

58+1 

5.60+0.03 

19/10/2013 

2456584.62 

19 000+400 

2.50+0.1 

2.0 

0.310+0.01 

320+30 

35+5 

58+1 

5.60+0.03 

20/10/2013 

2456585.67 

19 000+400 

2.47+0.1 

2.0 

0.320+0.01 

320+30 

50+5 

57+1 

5.59+0.03 

21/10/2013 

2456586.73 

19100+400 

2.50+0.1 

2.0 

0.280+0.01 

300+30 

60+5 

57+1 

5.60+0.03 

22/10/2013 

2456587.70 

19000+400 

2.40+0.1 

2.0 

0.285+0.01 

310+30 

60+5 

58+1 

5.60+0.03 

23/10/2013 

2456588.59 

19000+400 

2.40+0.1 

2.0 

0.285+0.01 

350+35 

60+5 

58+1 

5.57+0.03 

25/10/2013 

2456590.59 

18 650+400 

2.40+0.1 

2.0 

0.190+0.01 

300+30 

40+5 

58+1 

5.57+0.03 

26/10/2013 

2456591.59 

18600+400 

2.40+0.1 

2.0 

0.180+0.01 

270+30 

50+5 

57+1 

5.55+0.03 

27/10/2013 

2456592.59 

18570+400 

2.40+0.1 

2.0 

0.170+0.01 

270+30 

50+5 

56+1 

5.53+0.03 

28/10/2013 

2456593.60 

18570+400 

2.40+0.1 

2.0 

0.170+0.01 

270+30 

40+5 

56+1 

5.53+0.03 

30/10/2013 

2456595.71 

19000+400 

2.36+0.1 

2.2 

0.284+0.01 

330+30 

40+5 

57+1 

5.58+0.03 

31/10/2013 

2456596.71 

18600+400 

2.40+0.1 

2.0 

0.180+0.01 

270+30 

60+5 

57+1 

5.55+0.03 

01/11/2013 

2456597.57 

18600+400 

2.40+0.1 

2.0 

0.175+0.01 

260+30 

60+5 

57+1 

5.55+0.03 

07/11/2013 

2456603.59 

18570+400 

2.40+0.1 

2.0 

0.170+0.01 

270+30 

40+5 

56+1 

5.53+0.03 


(Petrov et al. 2014) As the mass-loss rate is a function of metal- 
licity and is derived from the strength of the Ha line, we expect 
a lack of uniqueness in wind models, which are based purely on 
the modeling of H transitions. 

This ambiguity is demonstrated in Fig. where we show 
in the top panel the Ha line computed for a constant mass-loss 
rate of 1.9 x lO^^Moyr^' and A(He)/A(H) ratios of 0.1, 0.2, 
and 0.4. To compensate for the vanishing Ha emission with in¬ 
creasing He abundance, we must increase the mass-loss rate to 
2.5 X 10“^ Mq yr“' and 3.5 x 10“^ Mq yr“', respectively. Then all 
three models produce approximately identical Ha profiles (bot¬ 
tom panel of Fig. [^. Consequently, uncertainties in the chemical 
composition (here He abundance) by a factor of may affect 
the derived mass-loss rate by a factor of 1.3-1.9. 

However, a higher He abundance also sensitively alters the 
strength of both the permitted and forbidden components of the 
He lines (see Fig. [^. When increasing the He abundance for a 
hxed mass-loss rate, the forbidden components develop promi¬ 
nent absorption features in the wings of the permitted lines. 
These manifest as pronounced asymmetries in the profiles of the 
permitted components, even if the projected rotational velocity 
is as high as ~ 60km s '. The presence and strength of the for¬ 
bidden component hence serve as useful diagnostics for the He 
surface abundance of slowly rotating stars. 

As the He line prohles observed in 55 Cyg do not show no¬ 
ticeable contributions from the forbidden components, a solar 
He abundance was found to be more appropriate to fit both the 
Ha and He lines. 

In relation to the modeling of the Si n and Si iii lines, we were 
able to obtain good fits using solar abundance, but the question 
remains of how sensitive the Si lines are to changes in their con¬ 
tent. The influence of non-solar abundance on was discussed 
in detail by [Lefever et al.| ( |2007[ ) and Markova & Puls (2008), 
who found that higher abundances result in lower values 
and vice versa. Such a change in temperature (ATeff > 1 000 K) 
would also affect the He lines. However, from our model fits, 
we do not find a discrepancy between the temperature obtained 
from the He and the Si lines. Therefore, the assumption of so¬ 
lar abundance for Si (within a range +0.15 dex) is well justified. 


All other elements are treated as background elements with solar 
abundance. 

4.2. Searching for pulsations 

As 55 Cyg resides inside the instability domain for BSGs (Saio 
et al. 2006) and a Cygni variables (Saio et al. 2013) , it is natu¬ 
ral to assume that pulsations might also be at work on the sur¬ 
face of this star. The strong variability in the Ha line (suggesting 
variable mass-loss) could indicate that the observed changes in 
the wind might be connected to pulsations detectable in photo- 
spheric line variations. To test this hypothesis, this section is de¬ 
voted to a detailed variability analysis of photospheric line pro¬ 
hles. 

4.2.1. Photospheric line profile variability 

In the absence of long-term photometric light curves, the 
presence of pulsations can also be tested by analyzing 
the variability of photospheric lines in high-quality opti¬ 
cal spectra. The spectral range observed with the Perek 
2m telescope contains several photospheric lines (see Fig.[^; 
the strongest is He 1/16678 A, followed by the deep ab¬ 
sorption lines of C11/1/16578,6583 A. Considerably weaker 
lines of Si 11/1/16347,6371 A, Ne i/l/l 6402,6507 A, and 
Nii/1/16482,6611 A are present as well. Despite the presence 
of these metal lines, we restricted most of our analysis to Her 
for two reasons. First, this line is not affected by telluric fea¬ 
tures, contrary to other lines such as Ha, C 11/1/16578,6583 A, 
Si II /l/l 6347,6371 A, and Nii /l/l 6482,6611 A. Any correction 
of the spectra with a telluric spectrum obtained from observing 
a standard star superimposes additional noise to the continuum 
level of the science spectra. Because He 1/16678 A falls outside 
regions of telluric pollution, no cleaning is required, and the ini¬ 
tial quality of the spectral line is preserved. Second, while only 
the C II lines are sufficiently strong for an analysis without suf¬ 
fering from an enhanced continuum noise level, these lines arise 
within a broad and shallow emission component centered on Ha 
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Table 2; As Tablej^ but for the observations taken in Ondfejov between 2009 and 2013. 


Date 

JD 

Teff 

log g 

P 

M 


f^macro(HQ') 

R* 

log LILq 

dd/mm/yr 


[K] 

[cgs] 


[10-^Moyr-i] 

[kms^'] 

[kms^'] 

[Re] 


15/08/2009 

2455059.44 

19000+1000 

2.40+0.1 

2.0 

0.290+0.01 

310+30 

35+5 

58+2 

5.60+0.03 

18/08/2009 

2455062.43 

19000+1000 

2.43+0.1 

2.0 

0.295+0.01 

340+35 

50+5 

54+2 

5.54+0.03 

23/08/2009 

2455067.49 

19000+1000 

2.43+0.1 

2.0 

0.340+0.01 

320+30 

30+5 

55+2 

5.55+0.03 

16/07/2010 

2455394.41 

18700+1000 

2.35+0.1 

2.0 

0.240+0.01 

250+25 

14+5 

63+2 

5.64+0.03 

31A/07/2010 

2455409.37 

18700+1000 

2.50+0.1 

2.0 

0.195+0.01 

250+25 

35+5 

57+2 

5.56+0.03 

31B/07/2010 

2455409.53 

18700+1000 

2.50+0.1 

2.0 

0.195+0.01 

250+25 

40+5 

57+2 

5.56+0.03 

OlA/08/2010 

2455410.35 

18700+1000 

2.50+0.1 

2.0 

0.195+0.02 

250+50 

80+10 

57+3 

5.56+0.06 

OlB/08/2010 

2455410.41 

18700+1000 

2.50+0.1 

2.0 

0.195+0.02 

250+50 

90+10 

57+3 

5.56+0.06 

OlC/08/2010 

2455410.45 

18700+1000 

2.50+0.1 

2.0 

0.195+0.02 

250+50 

100+10 

57+3 

5.56+0.06 

04 A/09/2010 

2455444.38 

18650+1000 

2.50+0.1 

2.0 

0.175+0.01 

270+30 

20+5 

56+2 

5.54+0.03 

04B/09/2010 

2455444.41 

18650+1000 

2.50+0.1 

2.0 

0.175+0.01 

270+30 

30+5 

56+2 

5.54+0.03 

04C/09/2010 

2455444.45 

18650+1000 

2.50+0.1 

2.0 

0.175+0.01 

270+30 

40+5 

56+2 

5.54+0.03 

05/09/2010 

2455445.31 

19000+1000 

2.41+0.1 

2.2 

0.280+0.01 

310+30 

50+5 

55+2 

5.55 +0.03 

06/09/2010 

2455446.52 

19000+1000 

2.41+0.1 

2.2 

0.280+0.01 

310+30 

60+5 

55+2 

5.55+0.03 

11/09/2010 

2455451.27 

18600+1000 

2.50+0.1 

2.0 

0.195+0.01 

250+25 

120+10 

57+2 

5.55+0.03 

21/09/2010 

2455461.34 

18600+1000 

2.40+0.1 

2.0 

0.180+0.02 

250+50 

30+5 

58+3 

5.56+0.06 

22/09/2010 

2455462.32 

18600+1000 

2.48+0.1 

2.0 

0.153+0.02 

235+50 

55+5 

58+3 

5.64+0.06 

23/09/2010 

2455463.38 

18700+1000 

2.42+0.1 

2.0 

0.236+0.01 

260+25 

50+5 

63+2 

5.56+0.03 

24/09/2010 

2455464.38 

19000+1000 

2.45+0.1 

2.0 

0.305+0.01 

330+35 

35+5 

54+2 

5.54+0.03 

06/03/2011 

2455627.64 

19000+1000 

2.43+0.1 

2.0 

0.335+0.01 

320+30 

30+5 

55+2 

5.55+0.03 

08/03/2011 

2455629.60 

19000+1000 

2.43+0.1 

2.0 

0.360+0.01 

320+30 

20+5 

55+2 

5.55+0.03 

19/03/2011 

2455640.63 

19000+1000 

2.43+0.1 

2.0 

0.370+0.01 

320+30 

100+10 

55+2 

5.57+0.03 

06/09/2011 

2455811.40 

18700+1000 

2.46+0.1 

2.0 

0.265+0.01 

300+30 

30+5 

58+2 

5.58+0.03 

25/09/2011 

2455830.41 

18900+1000 

2.45+0.1 

2.0 

0.400+0.01 

700+70 

40+5 

55+2 

5.55+0.03 

16/10/2011 

2455851.38 

19100+1000 

2.50+0.1 

2.0 

0.238+0.01 

350+35 

50+5 

53+2 

5.53+0.03 

26/02/2012 

2455984.63 

19000+1000 

2.50+0.1 

2.0 

0.460+0.01 

600+60 

30+5 

52+2 

5.50+0.03 

30/09/2012 

2456201.23 

18600+1000 

2.30+0.1 

2.0 

0.210+0.01 

230+25 

80+10 

65+2 

5.67+0.03 

03/10/2012 

2456204.26 

18600+1000 

2.30+0.1 

2.0 

0.200+0.02 

230+50 

30+5 

65+3 

5.67+0.06 

14/10/2012 

2456215.28 

19000+1000 

2.43+0.1 

2.0 

0.293+0.01 

340+35 

30+5 

54+2 

5.54+0.03 

14/11/2012 

2456246.24 

19100+1000 

2.46+0.1 

2.0 

0.200+0.02 

450+90 

70+5 

58+3 

5.61+0.06 

16/03/2013 

2456368.52 

19000+1000 

2.45+0.1 

2.0 

0.300+0.01 

300+30 

40+5 

52+2 

5.51+0.03 

21/04/2013 

2456404.47 

19000+1000 

2.45+0.1 

2.0 

0.305+0.01 

330+35 

40+5 

54+2 

5.54+0.03 

25/04/2013 

2456408.47 

19000+1000 

2.45+0.1 

2.0 

0.305+0.01 

330+35 

30+5 

54+2 

5.54+0.03 

19/07/2013 

2456493.37 

19000+1000 

2.40+0.1 

2.0 

0.295+0.01 

310+30 

35+5 

58+2 

5.60+0.03 

20/07/2013 

2456494.35 

18700+1000 

2.35+0.1 

2.0 

0.200+0.01 

250+25 

50+5 

63+2 

5.65+0.03 

21/07/2013 

2456495.33 

18700+1000 

2.35+0.1 

2.0 

0.200+0.01 

250+25 

65+5 

63+2 

5.65+0.03 

22/07/2013 

2456496.35 

18700+1000 

2.35+0.1 

2.0 

0.200+0.01 

250+25 

110+10 

63+2 

5.65+0.03 

23/07/2013 

2456497.52 

18700+1000 

2.35+0.1 

2.0 

0.200+0.01 

250+25 

80+10 

63+2 

5.65+0.03 

25/07/2013 

2456499.55 

18700+1000 

2.35+0.1 

2.0 

0.200+0.02 

250+50 

150+15 

63+3 

5.65+0.06 

30/07/2013 

2456504.35 

18700+1000 

2.35+0.1 

2.0 

0.200+0.01 

250+25 

55+5 

63+2 

5.65+0.03 

01/08/2013 

2456506.57 

18700+1000 

2.40+0.1 

1.6 

0.215+0.01 

250+25 

40+5 

60+2 

5.60+0.03 

02/08/2013 

2456507.34 

19000+1000 

2.43+0.1 

2.0 

0.295+0.01 

340+35 

80+10 

54+2 

5.54+0.03 

03/08/2013 

2456508.34 

19000+1000 

2.43+0.1 

2.0 

0.295+0.01 

340+35 

70+5 

54+2 

5.54+0.03 


Notes. Teff values are obtained from simultaneous fits to the H and He lines. The table lists only those models for which it was possible to derive 
the model parameters. Modeled and non-modeled Ha profiles are shown in Figs.|4|6| 


(see Fig.[T]). The shape of the line profiles of these C ii lines can 
thus be (possibly severely) altered by the normalization process, 
imprinting additional errors. These metal lines are, therefore, 
less favorable for the analysis of reliable line profile variabili¬ 
ties, but can be used as a secondary criterion to prove pulsation 
activity. 


4.2.2. Moment analyses 

For further analysis, we excluded all spectra that have exposure 
times longer than 1800 s. This leaves us with a total of 325 spec¬ 
tra, most with exposure times between 400 s and 900 s. 

To search for line profile variabilities in He 1/16678 A, we 
applied the moment method (Aerts et al. 1992; North & Pal- 
tani 1994) This method is ideal, as it allows distinguishing be¬ 
tween line profile variability caused by stellar spots or pulsa¬ 
tions. We computed the first three moments of He i /16678 A us¬ 
ing the prescription of Aerts et al. (1992) The first moment cor- 
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Fig. 7: Fits to the lines of Siii/l/l4128,4130 A, and 
Si III 44552 A in the Arizona spectra. Fits to the different line 
transitions are achieved using individual values of Wmic and timac 
(for details see Sect.|4. 1.3||. 


responds to the radial velocit}0 This parameter varies between 
+5kms~' and -21kms“' over the total observing period from 
August 2009 to October 2013 (see Fig.[T0|. The offset from a 


^ We follow the notation of Aerts et al. (2010a) in which the moments 
(not normalized and not corrected for systemic velocity) are denoted by 
nij, and the normalized, corrected moments by {v'). 


: H ALPHA 



Fig. 8: Top panel: Variation of the Ha line for different He abun¬ 
dances: A(He)/A(H) = 0.1 (solid line), A(He)/A(H) = 0.2 (dot¬ 
ted line), and A(He)/A(H) = 0.4 (dashed line) and the same 
mass-loss rate. Bottom panel: Similar Ha line profiles derived 
using the same A(He)/A(H) ratios, as in the previous plot, but 
with different mass-loss rates: 1.9 x lO^^Moyr^^ (solid line), 
2.5 X 10“^ Mq yr“' (dotted line), and 3.5 x 10“^ Mq yr“' (dashed 
line). 



Wavelength Wavelength (A) 


Fig. 9: Variation of the He i lines computed for a rotating star 
(y sin i = 60 km s“') with He contents of N{He)IN{H)- 0.1 
(solid line), N(He)IN(H)- 0.2 (dotted line) and N(He)IN(H)- 
0.4 (dashed line). The intensity of the He i forbidden components 
at A 4470 and A 4922 serves as an ideal indicator for the He abun¬ 
dance. 


symmetric variation around zero velocity is caused by the star’s 
intrinsic radial velocity, which we determine from this plot to be 
frad.star - -8kms“h We Corrected the data for this stellar mo¬ 
tion before we continued with our analysis. Hence, the observed 
maximum amplitude of the radial velocity variation in the line 
profile is about +13 km s“*. 

Next, we compared the behavior of the first and third mo¬ 
ments. Their time variability is shown in the middle and bottom 
panels of Fig. [m For better visibility and clarity, the plot is lim- 
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Fig. 10: Variation in the first moment (radial velocity, without 
correction for systemic velocity) of the He id 6678 A line in the 
period August 2009 to October 2013. 


Table 3: 55 Cyg periodicities 


No. 

/ 

(day 

dfx 10*" 
-') 

A dA 

(kms-') 

(h t/O 

(rad) 

/i 

0.044366 

0.03 

4.39 

0.03 

2.43 

0.01 

/2 

0.995570 

0.03 

3.31 

0.04 

3.71 

0.02 

h 

0.235716 

0.07 

2.54 

0.03 

3.44 

0.06 

/4 

0.974983 

0.17 

1.42 

0.03 

5.40 

0.09 

fs 

1.278593 

0.18 

1.37 

0.03 

3.35 

0.10 

fe 

2.400132 

0.17 

1.54 

0.02 

1.21 

0.11 

fi 

0.506733 

0.22 

1.30 

0.07 

5.55 

0.13 

h 

4.163120 

0.26 

0.87 

0.01 

5.78 

0.17 

/9 

6.994100 

0.21 

0.95 

0.01 

1.76 

0.14 

/lO 

3.104664 

0.18 

1.18 

0.02 

5.29 

0.10 

/ll 

0.188665 

0.42 

1.09 

0.04 

3.43 

0.26 

/l2 

1.462843 

0.40 

1.06 

0.03 

2.77 

0.25 

/l3 

0.296804 

0.24 

1.11 

0.03 

3.31 

0.14 

/l4 

0.652575 

1.21 

0.79 

0.03 

1.24 

0.77 

/l5 

1.181013 

1.55 

0.55 

0.04 

6.03 

1.00 

/l6 

0.617903 

0.30 

1.00 

0.02 

0.04 

0.17 

/l7 

2.593881 

0.45 

0.75 

0.02 

5.46 

0.28 

/l8 

9.003491 

0.52 

0.54 

0.01 

2.09 

0.34 

/l9 

3.780456 

1.10 

0.51 

0.02 

4.94 

0.79 


Notes. Phases are referenced to HID = 2455059.00 d. 


ited to a period of 25 days, ranging from 2013 July 19 to August 
12. This period contains time-series within seven nights, four at 
the beginning and three at the end. Obviously, the third moment, 
which describes the skewness (i.e., the asymmetry) of the line 
profile displays a strong time variability. In addition, the third 
moment varies in phase with the radial velocity, which is charac¬ 
teristic for stellar pulsations (see, e.g., Aerts et al. 1992) Similar 
behavior is seen in other seasons. 

Additional proof for pulsation activity in the atmosphere of 
55 Cyg is provided by comparing moments from different ele¬ 
ments. Although, as mentioned above, the other lines in the spec¬ 
trum have a lower quality, we computed the moments of the lines 
Cud 6578 A and Si ii d 6347 A, but restricted to the time-series 
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Fig. 11: Variation in EW (top), hrst (radial velocity, middle), and 
third moment (skewness, bottom) of He i d 6678 A in the observ¬ 
ing period 2013 July 19 to August 12. 


obtained during the nights 2013 July 21-22, which have the high¬ 
est S/N. The results for the hrst moment of these lines compared 
to He I d 6678 A are shown in Fig. 


12 


Obviously, all three lines 
vary identically. Such a behavior excludes the possibility of stel¬ 
lar spots. If the star has a patchy surface abundance pattern, as 
is common in chemically peculiar stars, for instance, the radial 
velocities of the He lines vary typically differently from those of 
the metal lines (e.g.. Briquet et al. 2001, 2004; Lehmann et al. 
2006) Instead, the identical variability of the line prohles of both 
He and metals in 55 Cyg is very typical for pulsations. 


The clear offset of the radial velocity curves of different lines 
could indicate that a velocity gradient exists within the atmo¬ 
sphere of 55 Cyg. In addition, the Si ii d 6347 A line, which is 
the weakest of these three, displays negative velocities, but the 
acceleration is in phase with the other two lines that display pos¬ 
itive velocitites. Unfortunately, the resolution of our spectra is 
too low to perform proper mode identihcation, so that no clear 
statements about the nature of the pulsations (radial, non-radial) 
can be made. 
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Fig. 12: Variation in the first moment (radial velocity) of 
He Id6678 A (dots). Cud6578 A (squares), and Si iid6347 A 
(triangles) during the nights 2013 July 21-22. 



0 2 _i6 8 10 

Frequency [d j 


Fig. 13: Lomb-Scargle periodogram of the original radial veloc¬ 
ity data (top) and of the data cleaned with the identified 19 sinu¬ 
soids (bottom). 

4.2.3. Frequency analysis 

The results from the moment analysis suggest pulsations as the 
most plausible cause for the observed line profile variability in 
He, as well as in the metal lines. The clear next step is to analyze 
the radial velocity curve to search for periods. We performed a 
frequency analysis using a similar strategy as the one described 
by Saio et al. (2006) We first computed a Lomb-Scargle peri¬ 
odogram (Scargle 1982) for frequency identification and ampli¬ 
tude estimation, which is displayed in the top panel of Fig.[T3| 
The periodogram finds no single dominant period. Instead, it 
seems that multiple significant peaks are present. Such a behav¬ 



1500 1510 

HJD-2455000 [d] 



Fig. 14: Fit (solid line) to the radial velocity variations (dots) of 
55 Cyg in the total period 2013 July 19 to August 12 (top) and 
zoomed in for the two marked series (bottom). 


ior is known for the a Cygni variables (Lucy 1976) for instance, 
and has also been reported for the periodograms of other BSGs 
(e.g., Saio et al. 2006) and we therefore used a fitting function 
consisting of a series of sinusoids. To refine the parameters for 
the identified frequency, we applied a nonlinear least-squares fit. 
A sine curve was fit to the data and subtracted. Then we searched 
for the next frequency within the residuals. In each subsequent 
cycle, the sum of all previously identified sine curves was sub¬ 
tracted from the original data. The procedure was stopped as 
soon as the S/N of the next frequency was too low. 

We define the noise level cr^es as the mean amplitude of the 
residual in a frequency interval around the identified frequency. 
The noise level in the frequency range 0 to 10 cycles day ' is 
0.13kms *. Frequencies with an amplitude A > 4crjes corre¬ 
spond to a confidence level of 99.9%, and are considered as clear 
detections. Those with A > 3.6crres correspond to a confidence 
level of 95%. In total, we identified 19 frequencies. Only the last 
one has an amplitude of 3.7crj.es, and should thus be considered 
as less accurate. All other detections have amplitudes A > 4cr^ss- 

With all frequencies identified, we again performed a least- 
squares fit using the parameters of all periods as first-guess val¬ 
ues. The final refined parameters of each sine curve with their er¬ 
rors are listed in Table]^ The Lomb-Scargle periodogram of the 
data, pre-whitened with the 19 identified periodicities, is shown 
in the bottom panel of Fig.[T3| 

Next, we computed the sum of the 19 individual sinusoids 
as a function of time (i.e., HJD) and compared it to the obser- 
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0 500 1000 1500 

HJD-2455000 [d] 


Fig. 15: Variation in EW (top) of the Hei/i6678A line over 
the total observing period, the radial velocity (middle) and the 
residuals of the radial velocity (bottom) after pre-whitening with 
the superposition of 19 sinusoids. 


vations. The synthetic radial velocity curve looks rather chaotic. 
For clarity and better visibility, we plot the synthetic curve lim¬ 
ited to a period of 25 days in the top panel of Fig.and over¬ 
plot the values obtained from the observations. In this plot, we 
mark two series, which correspond to the time interval where we 
obtained time series of observations of four and three consecu¬ 
tive nights. The comparison between the model and the observa¬ 
tions of these series is shown, zoomed in, in the bottom panels 
of Fig.[T4| 

In addition to the radial velocity of the He 1/16678 line, we 
also computed the line EW. Its value over the complete four- 
year observing period is displayed in the top panel of Fig.[^ It 
is found to be quite variable as well, which is an indication for 
temperature variations in the stellar atmosphere (e.g., De Ridder 
et al. 2002) We phased the data (both EW and first moment) to 
the longest period of 22.5 days identified from the radial velocity 
curve. This phase diagram is shown in Fig.[^ The radial veloc¬ 
ity clearly displays a sinusoidal behavior (bottom panel, with the 
sine curve of the identified period overplotted). We reveal vari¬ 
ability with the same period in EW (top panel), but with a phase 
shift of 94°. Considering that this period could belong to a ra¬ 
dial (strange) mode, a phase shift in EW of 90° would imply 
an adiabatic pulsation mode, while slight deviations are usually 
interpreted as non-adiabatic effects. Moreover, the large scatter 
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Fig. 16: Radial velocity and EW of the H1/16678 line phased to 
the longest period (22.5 d). The corresponding sine curve with 
the parameters from Table|^is included in the phase diagram of 
the radial velocity, and a tentative fit to the EW, which reveals 
the same period but is phase-shifted. 


in the phased EW can be explained by the multiperiodicity (De 
Ridder et al. 2002) 


5. Discussion 

We interpret the spectroscopic variability observed in the blue 
supergiant star 55Cyg in terms of photospheric oscillation 
modes and a time-variable radiation-driven wind. We discuss the 
reliability of the identified frequencies and the variation of wind 
parameter models. We also discuss the possibility that photo- 
spheric and wind variations might be related. 

5. 1. Reliability of the frequencies 

The agreement between the synthetic radial velocity curve and 
the observations seems to be quite good (see Fig.[l4|, but not 
all data fall on the model curve. This is also visible from the 
final residuals, which still display some scatter. These residuals 
are shown for the full observing period in the bottom panel of 
Fig.[T^ in comparison to the original observed radial velocity 
variation (middle panel). 

There are three possible reasons why we did not obtain better 
hts to the radial velocity curve. First, the number of observations 
is not sufficient to properly determine the stellar radial veloc¬ 
ity component. Hence, the data might not be exactly symmet¬ 
ric around zero, influencing the determination of the frequencies 
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and amplitudes of the periods. Second, not all data points have 
the same high quality, and there are a number of rather large gaps 
in the observations. Both lead to an increase in the noise level 
of the amplitude spectrum. Hence, from ground-based spec¬ 
troscopy we identify only 19 frequencies, while (many) more 
might be present. For comparison, 48 frequencies were iden¬ 
tified in the supergiant HD 163 899 by Saio et al. (2006) us¬ 
ing space-based photometry. Third, nothing is known about the 
stability of the periodic variations in 55 Cyg. As the observing 
period stretches over more than four years, some of the peri¬ 
odicities might have changed or completely disappeared, while 
others might have appeared. Such a behavior of flipping peri¬ 
odicities was reported by Markova et al. (2008) for the BSG star 
HD 199478. Furthermore, Aerts et al. (2010b) recognized a sud¬ 
den change in the amplitude of the pulsation mode in the BSG 
star HD 50 064. 

Consequently, precise determinations of currently present 
periodicities in 55 Cyg must be obtained from a better dataset, 
ideally from continuous spectroscopic observations over a long 
time interval, combined with a continuous photometric light 
curve. This is a non-trivial task, however. 

5.2. Stellar rotation 

The presence of multiple periodicities with significant velocity 
amplitudes in the photospheric lines implies that their profiles 
are affected in two ways. First, the lines can be expected to be 
broadened (referred to as macroturbulence) in excess of the typ¬ 
ical broadening mechanisms acting in stellar atmospheres, in¬ 
cluding stellar rotation. And second, the profile shape is usually 
not symmetric anymore. The latter is confirmed from our mo¬ 
ment analysis (see Fig. [n]). However, the amounts of both rota¬ 
tion and macroturbulence are not known initially. The determi¬ 
nation of proper stellar rotation rates is particularly important be¬ 
cause rotation plays a major role in the evolution and mass-loss 
behavior of massive stars (see, e.g., Maeder & Meynet 2000; 
Meynet & Maeder 2005; Ekstrom et al. 2012; Simon-Diaz & 
Herrero 2014) 

Investigations based on different methods, such as the 
Fourier transformation (FT) or the goodness-of-fit (GOF), are 
usually applied with the aim to determine proper u sin i values 
and the contribution of macroturbulence. However, a degener¬ 
acy exists, which allows fitting observed line profiles with a di¬ 
versity of different combinations of vsini and Wmacro (e-g-, Ryans 
et al. 2002) This degeneracy was also found for 55 Cyg by Ju- 
rkic et al. (2011) from modeling silicon lines in lUE spectra. 
Eurthermore, the major requirement for a reliable result from the 
ET method is a symmetric signal. Application to asymmetric line 
profiles, as is typically expected for pulsating stars, immediately 
delivers false results. Studies of the effect of pulsational activ¬ 
ity in stellar atmospheres on resulting line profiles revealed that 
the collective effect of multiple periodicities, for example, in the 
form of low-amplitude non-radial gravity modes, can cause sub¬ 
stantial line broadening, which might be identified with macro¬ 
turbulence (Aerts et al. 2009; Simon-Dfaz et al. 2010a) These 
investigations also revealed that the ET method fails as soon as 
pulsations play a non-negligible role in the line width. More¬ 
over, Aerts et al. (2014) showed that results from the ET of line 
profiles of known pulsators vary substantially over the pulsation 
cycle. 

Even though our observations have extremely good S/N, 
which is considered as a major prerequisite for reliable results 
from the ET (Simon-Dfaz & Herrero 2007), it is not possible 
to determine consistent v sin i values for 55 Cyg. This is shown 



40 60 80 40 60 80 

V sin i (km s ) v sin i (km s ) 


Eig. 17; Top; Selection of values of ysinf obtained with the ET 
method for He 1/16678 A (left) and Cii/16578 A (right) within 
the night of 2013 July 22. Only results with a clear first minimum 
are shown. Bottom; Range in v sin i values obtained within three 
consecutive nights. 


in Fig.[T7j where we plot the results from the ET analysis of 
He I /16678 A and C ii /16578 A during different times within one 
night (top panel) and extreme values obtained during three con¬ 
secutive nights (bottom panel). Eor many line profiles (for both 
He I and C n) the ET of our data deliver no clear minimum and, 
hence, no value for v sin i could be found. 
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From the profiles with clear zero points in their FT, C ii sug¬ 
gests usini values between 44kms“' and 48kms“' (which re¬ 
semble the rotational values of 40 km s“^ and 45 km s“^ obtained 
for the Sin lines). The values obtained from the Hei line are 
in most cases significantly higher and with a higher spread in 
values, for example, from about 40km s“* to about 65kms“'. 
While a shift in the position of the first zero point and the appear¬ 
ance of saddle points and local minima in the FT can be expected 
in the case of pulsations (Aerts et al. 2014), the discrepancy in 
the values obtained from the two investigated lines is difficult 
to understand as the FT method usually results in very similar 
II sin i values when applied to He and metal lines (Simon-Dfaz & 
Herrero 2007) While the variation in values obtained from C ii 
is small, it may still imply that macroturbulent line broadening 
(and hence the underlying pulsations) influences the position of 
the zero points. In that context, strongly varying vsini values 
within a single night were obtained from the FT method applied 
to the Si 11/16347 line of the BSG HD 202 850 by Tomic et al. 
(2015), 


In addition to the seemingly higher usin/ value resulting 
from the FT, the He 1/16678 A line also appears to be special 
in other aspects. For instance, the measured radial velocity is 
higher than in the Cii/i6578A line (see Fig. 12 1 . Moreover, 
He 1 46678 A is the only line in our sample (apart from Ha) 
that requires substantially higher values for both the micro- and 
macroturbulent velocities than all other modeled photospheric 
lines. A possible explanation could be that the Cii/16578 A and 
the He 1/16678 A lines are formed in different regions in the at¬ 
mosphere with Cii originating in deeper layers, in agreement 
with the Cii line being much weaker than the Hei line (see 
Fig.[^. The line-forming region for He i might be substantially 
larger and may be extending into the base of the wind, where 
higher turbulent velocities (i.e., higher than in the photosphere) 
might be expected due to instabilities, possibly resulting from 
coupling of the pulsation modes. 

Although the FT analysis of the line profiles has severe prob¬ 
lems in the presence of strong macro turbulence and possible 
wind contamination of the profiles, the range of values obtained 
from both lines is on the same order as the value of 55 + 5 km s ' 
found from the line profile fits (see Sect. |4.1.3] l. However, as long 
as the pulsation modes cannot be identified, the individual con¬ 
tributions from the two different, but physically very important 
mechanisms acting on the line profiles, that is, rotation and pul¬ 
sations, remain uncertain. 


5.3. Connection between pulsations and mass-loss 

Our analysis revealed multiple periods, which can be interpreted 
as stellar pulsations. Hence, the question arises whether these 
pulsations might be suitable to trigger time-dependent enhanced 
mass loss, a question that was also addressed by Lefever et al. 
(2007) for a large sample of BSGs with no firm answer. 

Two frequencies in the set of frequencies identified in 
Sect. |4.2.3] f 3 (4.2 days) and fn (5.3 days), are close to the spec¬ 
troscopic period of 4-5 days reported by Granes (1975) and the 
periodicity of 4.88 days found in the photometric light curve by 
Koen & Eyer (2002) Twelve frequencies correspond to oscil¬ 
lations with periods spreading from ~ 2.7 hours to ~24 hours. 
The frequencies fs and fe might be harmonics of fis and fie, 
respectively. Meanwhile, the frequencies associated with peri¬ 
ods of about one day (f 2 and f 4 ) might be related to the daily 
variations observed in the intensity of the Ha line because their 
similar frequencies cause beats. 
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We would like to stress that the daily changes in line profile 
shape and intensity seen in our data should not be interpreted as 
a real change in mass loss. Each of our spectra has been mod¬ 
eled independently, considering a spherical smooth wind with a 
yS velocity law. However, any density enhancement at the base 
of the wind due to some mass ejection (triggered by pulsation 
or any other instability) will need time to travel through the Ha 
line-forming region and hence contribute to the Ha line profile 
in more than just one night. The night-to-night variations that we 
see in the data can therefore not be regarded as real changes in 
the wind parameters (mass loss and terminal velocity) but should 
rather be interpreted as rapid changes in the wind density, that is, 
wind inhomogeneities, traveling through the wind and influenc¬ 
ing the Ha line profile. 

It is also worth mentioning that there were three occasions 
at which the Ha line completely (or almost completely) dis¬ 
appeared from the spectrum, indicating that in these particular 
cases the wind emission was reduced so that it compensated by 
coincidence for the photospheric absorption line. The first time 
noted was on 2010 July 2-4 by Maharramov (2013), and in our 
own data we recognized it on 2012 July 24 and in 2013, around 
July 26. Interestingly, after the re-appearance of the Ha in emis¬ 
sion, its profile does not recover the shape seen before the can¬ 
cellation. Instead, a strong, pure emission evolves during 2013 
August 10-12, which might be explained as due to an enhanced 
amount of material appearing in the Ha line-forming region. 
Eurthermore, inspection of the behavior of the He 1/16678 line 
in the same period reveals a strong increase in the radial velocity 
(middle panel of Fig. 11 1 with a simultaneous drop in EW (top 
panel of Eig. M- 

If the occasions of annihilation of emission and absorption 
components of the Ha line were periodic, the possible period 
would be on the order of about one year. We may speculate that 
the wind conditions of 55Cyg may have changed some years 
ago, because no earlier observational work on 55 Cyg mentions 
phases in which the wind emission replenishes the photospheric 
absorption line. 

On the other hand, according to the stellar parameters de¬ 
rived in this work, log L/Lq = 5.57, and the current stellar 
mass of 34 Mq, we find a ratio for LjM ~ 1.1 x 10"'Lo/Mo , 
which is high enough for the occurrence of radial strange modes 
(Gautschy & Glatzel 1990; Glatzel 1994; Saio et al. 1998) The 
stellar properties of 55 Cyg resemble those of HD 50 064, for 
which Aerts et al. (2010b) argued that the star is building up 
a circumstellar envelope, as seen in the variable Balmer lines, 
while pulsating in a radial strange mode. A similar scenario 
might also hold for 55 Cyg. The first frequency, fi (22.5 days), 
has the proper length to be interpreted as radial strange mode 
as predicted from theory (Saio et al. 2013) This period might 
tentatively be connected to the mid-term variation observed in 
the Ha line and the mass-loss rates. Considering the sequence of 
observations obtained in October 2013 in Arizona, it seems that 
the mass-loss rate may keep more or less a constant value over a 
time interval of roughly a week, followed by a jump to a lower 
mass-loss regime. This second regime might have lasted seven or 
perhaps more days (neglecting the sudden jump in the mass-loss 
observed on October, 30). Another period of at least five days of 
constant behavior can be recognized in July 2013. Whether this 
(quasi-)periodic variability in mass-loss might be recovered in 
other seasons cannot be asserted because of the often large gaps 
in observations. 

Additional support for a possible connection between pul¬ 
sations and enhanced mass loss is provided by the variable line 
shapes of Ha and by the faint shallow symmetric wings extend- 
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ing to +1200 km s * (described in Sect.[^. The velocity of these 
wings is definitely higher than the wind terminal velocities de¬ 
termined for 55 Cyg, so that they are not formed by the current 
wind material. Instead, we interpret these wings as due to in¬ 
coherent electron scattering in the outer layers of the star’s ex¬ 
tended envelope above the line formation region. As the wings 
were not reported by other observers (with one exception, see 
Ebbets 1982), but are rather stable during the past four years, the 
pulsational activity leading to the formation of an optically thick, 
expanding atmosphere might have started or increased within the 
past few years. 

5.4. Evolutionary state of 55 Cyg 

According to the derived effective temperature and luminosity, 
we might connect the star with an evolutionary track of 40 Mq. 
The assignment of a proper evolutionary state of the star is ham¬ 
pered by the fact that massive stars can cross the temperature 
range of BSGs more than once during their evolution. The loca¬ 
tion of a star in the HRD alone is thus not sufficient, and addi¬ 
tional criteria are needed for an unambiguous evolutionary state 
determination. One such criterion is provided by the recent cal¬ 
culations of the occurrence of possible pulsation modes during 
the post-main sequence evolution of massive stars by Saio et al. 
(2013) Their investigations revealed that during the post-RSG 
evolution, BSGs can pulsate in many more modes and, in partic¬ 
ular, in radial strange modes, which is not possible while the star 
is on its first redward evolution. Moreover, the observations of 
pulsation properties with very many excited modes are possible 
only after the star has experienced a strong mass-loss during the 
RSG phase. 

The numerous periodicities identified in 55 Cyg display a 
mixture of p- and g-modes. According to computations of pul¬ 
sation modes in BSGs, for rotating stellar models in the range of 
initial masses of 20-25 Mgyr^' (Saio et al. 2006, 2013), the co¬ 
existence of both types of modes in the pre-RSG stage is found 
for stars with temperatures T^g ^ 20 000 K, while cooler stars 
are expected to display very few pulsations. Even though the 
temperature of 55 Cyg, T^g - 19 000 K, is slightly lower, its pe¬ 
riodicities are too numerous to fit in a pre-RSG scenario. Con¬ 
sequently, a post-RSG stage is more plausible mainly because 
of the multiple pulsation modes (with periods 30 days) that 
we find in this star. Eurthermore, the classification of 55 Cyg as 
a post-RSG is supported by the presence of the bow shock fea¬ 
ture detected by the Wide-held Infrared Survey Explorer (WISE, 
Wright et al. 2010) in band 4 at 24 /rm (see Eig. [TS] ), which 
formed in the past by the interaction of a strong wind phase with 
the ISM. 

But this post-RSG scenario does not seem to ht the low He 
surface abundance we derived from our line prohle modeling. 
This apparent contradiction was also found in at least two other 
BSGs (Rigel and Deneb, see Saio et al. 2013) To solve this prob¬ 
lem, Georgy et al. (2014) recently showed that new stellar evolu¬ 
tion calculations, using the Ledoux criterion for convection, rec¬ 
oncile the expected pulsational properties of a post-RSG phase 
with lower surface He chemical abundance (values that were pre¬ 
viously expected only during the pre-RSG phase assuming the 
Schwarzschild criterion for convection). 


6. Conclusions 

55 Cyg shows both long- and short-term variability in its spectro¬ 
scopic and photometric data. WISE IR images reveal a large bow 



Eig. 18: Bow shock features detected by WISE in band 4 at 24 
yum. The size of the image is 10' x 10'. 


shock structure originating from the interaction of a strong mass- 
loss event with the ISM. This strong wind phase is expected to 
occur during evolution through the RSG stage. 

To understand the origin of the variability observed in the 
photospheric lines, we performed a moment analysis based on 
almost 400 optical spectra spread over about five years and iden¬ 
tified 19 periods (although many more might be present). Twelve 
of these periods spread from ~ 2.7 hours to ~ 24 hours and might 
belong to p-modes, and 6 are between 1 and 6 days long and 
could represent g-modes. The longest period of 22.5 days could 
be related to the theoretically predicted strange mode pulsations. 
Support for this interpretation comes from the value of the L/M 
ratio, which is of the correct order of magnitude for strange mode 
excitation. The existence of many short pulsation modes might 
influence the line profile of the photospheric lines, and might, 
hence, be responsible for the large uncertainties in the derived 
rotational velocity of the star and the variety of macroturbulent 
velocities measured in different lines. 

We modeled the H, He i. Si ii, and Si iii lines using the NLTE 
atmosphere code EASTWIND and derived stellar and wind pa¬ 
rameters. We found that the star undergoes episodes in which the 
determined mass-loss can vary by a factor of 1.7-2 on timescales 
of about two to three weeks. We also detected a long-term period 
in which the star exhibits no H emission lines. In the Ha line pro¬ 
file, we detected strong night-to-night intensity and profile shape 
variability, indicating rapid changes in the wind that might be in¬ 
terpreted as local material enhancements, caused by short-term 
mass ejection events traveling through the Ha line-forming re¬ 
gion. Multiple periods in both the wind and photospheric ma¬ 
terial might suggest a possible link between pulsational activ¬ 
ity and enhanced mass-loss in 55 Cyg. Moreover, combining the 
high derived spectroscopic mass of M* = 34 Mq, the many ex¬ 
cited pulsation modes, the evidence of a huge past mass-loss 
phase, and the presence of CNO processed material at the sur¬ 
face, 55 Cyg is most probably in a post-RSG phase, crossing the 
HRD toward the blue region. As such, the star conforms to char¬ 
acteristics resembling classical aCygni variables (Georgy et al. 
2014; Saio et al. 2013) 
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